Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 4 February 2008 (MN KTpX style file v2.2) 



Oh 

00 



Reconstructed Density and Velocity Fields from the 2MASS 
Redshift Survey 

Pirin Erdogdu 1 ' 2 ' 3 , Ofer Lahav 3 ' 4 , John P. Huchra 5 *, Matthew Colless 6 , Roc M. Cutri 7 , 
Emilio Falco 5 * Teddy George 6 , Thomas Jarrett 8 , D. Heath Jones 7 , Lucas M. Macri 9 , 
Jeff Mader 10 , Nathalie Martimbeau 5 , Michael A. Pahre 5 , Quentin A. Parker 711 , Anais 
^ ■ Rassat 4 , Will Saunders 7 

C ~ 3 1 Department of Physics, Middle East Technical University, 06531, Ankara, Turkey 

£Nj 2 School of Physics & Astronomy, University of Nottingham, University Park, Nottingham, NG7 2RD, UK 

3 Institute of Astronomy, Madingley Road, Cambridge, CB3 0HA, UK 

4 Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK 

5 Harvard-Smithsonian Centre for Astrophysics, 60 Garden Street, MS-20, Cambridge, MA 02138, USA 

6 Canada-France-Hawaii Telescope Corporation, 65-1238 Mamalahoa Hwy Kamuela, Hawaii 96743, USA 
' 7 Anglo-Australian Observatory, PO Box 296, Epping,NSW 2052, Australia 

£S1 ' 8 Infrared Processing and Analysis Center, California Institute of Technology, Pasadena, CA 91125, USA 

9 National Optical Astronomy Observatory, 950 North Chrry Avenue, Tucson, AZ 85726, USA 
7-H ■ 10 W.M. Keck Observatory, Kamuela, HI 96743, USA 

^ ' 11 Department of Physics, Macquarie University, Sydney, NWS 2109, Australia 
I/~) ' * Guest Observer, Cerro Tololo Interamerican Observatory, operated by AURA for the National Science Foundation 

o" 
o 
o . 

1—1 , 4 February 2008 

O ; 

J3 '. ABSTRACT 

We present the reconstructed real-space density and the predicted velocity fields from the 
1. \ Two Mass Redshift Survey (2MRS). The 2MRS is the densest all-sky redshift survey to 

date and includes about 23,200 galaxies with extinction corrected magnitudes brighter than 
' K s = 11.25. Our method is based on the expansion of these fields in Fourier-Bessel fianc- 

ee , tions. Within this framework, the linear redshift distortions only affect the density field in 

L* ■ the radial direction and can easily be deconvolved using a distortion matrix. Moreover, in 

. 5^ | this coordinate system, the velocity field is related to the density field by a simple linear 

, transformation. The shot noise errors in the reconstructions are suppressed by means of a 

5_i ■ Wiener filter which yields a minimum variance estimate of the density and velocity fields. 

Using the reconstructed real-space density fields, we identify all major superclusters and 
voids. At 50 h^ 1 Mpc, our reconstructed velocity field indicates a back-side infall to the 
Great Attractor region of Ui n f a u = (491 ± 200)(/3/0.5) kms -1 in the Local Group frame 
and Ujnfall = (64 ± 205)(/3/0.5) kms -1 in the cosmic microwave background (CMB) frame 
and (3 is the redshift distortion parameter. The direction of the reconstructed dipole agrees 
well with the dipole derived by Erdogdu et al. (2006). The misalignment between the recon- 
structed 2MRS and the CMB dipoles drops to 13° at around 5000 kms -1 but then increases 
at larger distances. 

Key words: astronomical data bases:surveys-cosmology: observations - large-scale structure 
of universe - galaxies: distances and redshifts 



1 INTRODUCTION 

A cosmographical description of the galaxy distribution is crucial 
to our understanding of the mechanisms of structure formation that 
generate the complex pattern of sheets and filaments comprising 
the 'cosmic web'. Today, there are many more redshifts available 
for galaxies than velocity measurements. This discrepancy has in- 



spired a great deal of work on methods for the reconstruction of 
the mass and velocity fields from redshifts alone. These methods 
use a variety of functional representations (e.g. Cartesian, Fourier, 
spherical harmonics or wavelets) and smoothing techniques (e.g. 
a Gaussian or a top-hat sphere). However, several complications 
arise when one tries to reconstruct peculiar velocity and real-space 
mass density fields directly from redshift surveys. First, one must 
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relate the light distribution to the mass field. On large scales, this 
relationship is generally assumed to be linear with a proportionality 
constant (but see e.g. Wild et al. 2005; Conway et al. 2005; Mari- 
noni et al. 2005 who argue that biasing is non-linear and redshift 
dependent even on large scales). Another common assumption is 
that the galaxy distribution samples the underlying smooth mass 
field and the finite sampling of the smooth field introduces Pois- 
son 'shot noise' errors 1 . Yet another important physical problem 
associated with the recovery of the real-space mass density field 
from the galaxy density field is the correction of the distortions in 
the redshift space clustering pattern. Furthermore, in redshift sur- 
veys the actual number of galaxies in a given volume is larger than 
the number observed, particularly in flux limited samples where at 
large distances only the very luminous galaxies can be seen. 

In this paper, we recover the real-space density and predicted 
velocity fields of the 2MRS (Huchra et al. 2005). The reconstruc- 
tion procedure used is based on linear theory and closely follows 
Fisher et al. (1995, hereafter FLHLZ). The density field in redshift- 
space is decomposed into spherical harmonics and Bessel functions 
(Fourier-Bessel functions). Then the real-space density and veloc- 
ity fields are reconstructed from those in the redshift- space using a 
Wiener filter. 

The idea of cosmography using spherical harmonics goes back 
to Peebles (1973) but the method only gained popularity with the 
advent of all-sky surveys. Since then, it has been applied to the 
Infrared Astronomical Satellite (IRAS) selected galaxy surveys 
(e.g. to the IRAS 1.2 Jy Survey by Fabbri & Natale 1990; Scharf 
et al. 1992; Scharf & Lahav 1993; Fisher, Scharf & Lahav 1994; 
Lahav et al. 1994; Nusser & Davis 1994; Heavens & Taylor 1995 
and to IRAS Point Source Catalogue Redshift Survey (PSCz) by 
Tadros et al. 1999, Hamilton, Tegmark & Padmanabhan 2000; Tay- 
lor et al. 2001) and the peculiar velocity catalogues (e.g. to elliptical 
galaxy sample by Regos & Szalay 1989; to Mark III galaxy sam- 
ple by Davis, Nusser & Willick 1996 and Hoffman et al. 2001). 
The spherical harmonics technique has proven to be very conve- 
nient for addressing many of the problems inherent to the analyses 
of redshift surveys. The spherical harmonics are one of the most 
convenient ways of smoothing the noisy data. A careful choice 
of spherical coefficients greatly reduces the statistical uncertain- 
ties and errors introduced by the non-linear effects. Furthermore, 
the treatment of linear redshift distortions are particularly straight- 
forward in Fourier-Bessel space, especially for surveys which have 
almost all-sky coverage and dense sampling such as 2MRS. 

The Wiener-filtered fields are optimal reconstructions in the 
sense that the variance between the derived reconstruction and the 
underlying true field is minimised. As opposed to ad hoc smooth- 
ing schemes, the smoothing due to the Wiener filter is determined 
by the data. In the limit of high signal-to-noise, the Wiener filter 
modifies the observed data only weakly, whereas it strongly sup- 
presses the contribution of the data contaminated by shot noise. The 
three-dimensional Wiener reconstruction of cosmological fields in 
Fourier-Bessel space was first developed by FLHLZ and applied 
to the IRAS 1.2 Jy survey. Webster et al. (1997, hereafter WLF) 
expanded on the initial results of FLHLZ. Later, Schmoldt et al. 



(1999) applied the same method to the IRAS PSCz survey to re- 
cover density and velocity fields. 

This paper is structured as follows: in the next section, we de- 
scribe the Two Micron Redshift Survey. In Section 3, we give a 
brief summary of the decomposition of the density field in redshift- 
space into a set of orthogonal spherical harmonics and Bessel func- 
tions. In Section 4, we discuss how these redshift-space harmonics 
are related to those in real-space. In Section 5, we give expres- 
sions for the Wiener filter in Fourier-Bessel space. The implemen- 
tation of the method to 2MRS is given in Section 6, along with 
the maps of the density and velocity fields. In the same section, we 
also show reconstruction of the Local Group (LG) acceleration de- 
termined directly from harmonic coefficients. The last section con- 
tains a discussion and summary of the work. Throughout this paper, 
we assume that the galaxy distribution is related to the underlying 
smooth mass density field by the linear bias parameter, b, so that 
5 g (x) = b& m (x) . With this notation, the constant of proportional- 
ity between velocity and the density field is (3, where j3 — Q m os /b. 
It is worth emphasising that one of the ultimate goals of this work 
is to test linear biasing by comparing the reconstructed velocity 
fields with the observed peculiar velocity field. The distances are 
given in units of velocity and the Hubble constant is Ho = 100 h 
kms _1 Mpc _1 . 



2 THE TWO MICRON ALL-SKY REDSHIFT SURVEY 

The Two Micron All-Sky Redshift Survey (2MRS) is the densest 
all-sky redshift survey to date. The galaxies in the northern celestial 
hemisphere are being observed by the Fred Lawrence Whipple Ob- 
servatory (FLWO) 1.5-m telescope, the Arecibo 305-m telescope 
and the Green Bank 100-m Telescope. In the southern hemisphere, 
most galaxies at high galactic latitude (about 6000 galaxies) were 
observed as a part of the 6-degree field galaxy survey 2 (6dFGS, 
Jones et al. 2004; Jones et al. 2005) and the low galactic latitude 
galaxies are being observed at CTIO (by L. Macri and J.P. Huchra). 
The first phase of 2MRS is now complete. In this phase we obtained 
redshifts for approximately 23,150 2MASS galaxies from a total 
sample of 24,773 galaxies with extinction corrected magnitudes 
(Schlegel, Finkbeiner & Davis 1998) brighter than K s = 11.25. 
This magnitude limit corresponds to a median redshift of z ~ 0.02 
(« 6000 km s -1 ). Almost all of the 1600 galaxies that remain 
without redshifts are at very low galactic latitudes (|6| < 5°) or 
obscured/confused by the dust and the high stellar density towards 
the Galactic Centre. Figure 1 shows all the objects in 2MRS in 
Galactic Aitoff Projection. Galaxies with z ^ 0.01 are plotted in 
red, 0.01 < z < 0.025 are plotted in blue, 0.025 < z < 0.05 
are plotted in green and z ^ 0.05 are plotted in magenta. Galax- 
ies without measured redshifts are plotted in black. 2MRS can be 
compared with the deeper 2MASS galaxy catalogue (K< 14th mag) 
shown in Jarrett (2004, Figure 1). 



1 This is only an approximation. A more realistic model for galaxy clus- 
tering is the halo model (e.g. Cooray & Sheth 2002) where the linear bias 
parameter depends on the mass of the dark matter haloes where the galaxies 
reside. For this model, the mean number of galaxy pairs in a given halo is 
usually lower than the Poisson expectation. 



2 6dFGS has a much fainter magnitude limit (K s = 12.75) than 
2MRS and contains around 150,000 galaxy redshifts. The 6dFGS 
documentation and the database can be found on the WWW at 
http://www.aao.gov.au/local/www/6df/ 
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Figure 1. All Objects (24 788 galaxies) in the 2MASS Redshift Catalogue in Galactic Aitoff Projection. The plot contains Galaxies with z sC 0.01 plotted in 
red, 0.01 < z ^ 0.025 are plotted in blue, 0.025 < z < 0.05 are plotted in green and z 0.05 are plotted in magenta. Galaxies without measured redshifts 
are plotted in black. The zone of avoidance is outlined by dashed lines. 



2.1 Survey Completeness 

2MASS 3 has good photometric uniformity and an unprecedented 
integral sky coverage. The photometric uniformity is better than 4% 
over the sky including the celestial poles (e.g. Jarrett et al. 2000a, 
2003) and at magnitudes brighter than K s = 12, 2MASS is es- 
sentially complete down to very low galactic latitudes (to b 2 o, e.g. 
Huchra et al. 2005). The uniform completeness of the galaxy sam- 
ple is slightly limited by the presence of foreground stars, for a 
typical high latitude sky less than 2% of the area is masked. These 
missing regions are accounted for using a coverage map, defined as 
the fraction of the area of an 8' x 8'pixel that is not obscured by stars 
brighter than 10th mag. Galaxies are then weighted by the inverse 
of the completeness, although the analysis is almost unaffected by 
this process as the completeness ratio is very close to one for most 
parts of the sky. 



The stellar contamination of the catalogue is low and is re- 
duced further by manually inspecting the objects below a helio- 
centric redshift of cz = 200 km s" 1 . The foreground stellar con- 
fusion is highest at low Galactic latitudes, resulting in decreasing 
overall completeness of the 2MASS catalogue (e.g. Jarrett et al. 
20006) and consequently the 2MRS sample. Stellar confusion also 
produces colour bias in the 2MASS galaxy photometry (Cambresy, 
Jarrett & Beichman 2005) but this bias should not be significant for 
2MRS due to its relatively bright magnitude limit. 

2.2 Magnitude and Flux Conversions 

2MRS uses the 2MASS magnitude K20, which is defined 4 as the 
magnitude inside the circular isophote corresponding to a surface 
brightness of fiK s = 20mag arcsec -2 (e.g. Jarrett et al. 2000a). 



3 2MASS database and the full documentation are available on the WWW 
at http//w w w. ipac .caltech.edu/2mass. 
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The isophotal magnitudes underestimate the total luminosity by 
10% for the early-type and 20% for the late-type galaxies (Jar- 
rett et al. 2003). Following Kochanek et al. (2001, Appendix), an 
offset of A = —0.20 ± 0.04 is added to the K20 magnitudes. 
The galaxy magnitudes are corrected for Galactic extinction us- 
ing the dust maps of Schlegel, Finkbeiner & Davis (1998) and an 
extinction correction coefficient of Rk = 0.35 (Cardelli, Clay- 
ton & Mathis 1989). As expected, the extinction corrections are 
small for the 2MRS sample. The K s band ^-correction is derived 
by Kochanek et al. (2001) based on the stellar population models 
of Worthey (1994). The k-correction of k(z) = -6.0 log (1 + z), 
is independent of galaxy type and valid for z < 0.25. 

The flux, S, for each galaxy is computed from the apparent 
magnitudes using 

S = S(0mag)10- a4(K20+ZPO) , (1) 

where the zero point offset is ZPO = 0.017 ± 0.005 and 
5(0mag) = 1.122 x 10~ 14 ± 1.891 x 10* 16 Wcm~ 2 for the 
K B band (Cohen, Wheaton & Megeath 2003). 
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Figure 2. Redshift histogram for 2MRS galaxies and a least squares fit 
(Equation 2) to the data (black). For comparison, also plotted is a redshift 
histogram for PSCz galaxies (Saunders et al. 2000) (red). 



2.3 The Redshift Distribution and the Selection Function 



The redshift distribution of 2MRS is shown in Figure 2. The IRAS 
PSCz survey redshift distribution (Saunders et al. 2000) is also plot- 
ted for comparison. 2MRS samples the galaxy distribution better 
than the PSCz survey out to cz — 15000 kms" 1 . The selection 
function of the survey (i.e. the probability of detecting a galaxy as 
a function of distance) is modelled using a parametrised fit to the 
redshift distribution: 



dN(z) = Az 1 exp 



dz 



(2) 



with best-fit parameters of A = 116000 ± 4000, a = 2.108 ± 
0.003, 7 = 1.125±0.025 and z c = 0.025±0.001. This best-fit is 
also shown in Figure 2 (solid line). The overall selection function 
4>s{r) is the redshift distribution divided by the volume element 



Mr) 



dN 

dz 



dz\ 
dr), 



(3) 



where f2 s (~ 4-7T steradians) is the solid angle of the survey and r is 
the comoving distance. 



3 EXPANSION OF THE DENSITY, VELOCITY AND 
GRAVITATIONAL POTENTIAL FIELDS IN 
SPHERICAL HARMONICS 

This section describes the spherical harmonic expansion of the den- 
sity and velocity fields. Spherical harmonics are commonly em- 
ployed in cosmography of whole-sky surveys for two main rea- 
sons. Firstly, separating the observed density field into angular and 
radial modes concentrates the linear scale redshift distortions into 
one dimension, making distortions easier to deconvolve. Secondly, 
both spherical harmonic and Bessel functions are orthogonal and 
together they form eigenfunctions of the Laplacian operator. These 
properties lead to very simple relationships between the density, 
the velocity and the potential fields in harmonic space, allowing the 
straightforward reconstruction of one field from another observed 
field. The method used in this paper is based on that developed by 
FLHLZ. Therefore, we do not give a full description of our method, 
we only outline the formalism. 



3.1 Expansion of the density field 

The density field, p(r), can be expanded as a Fourier-Bessel series 
in the following way: 

! m „ +; n ma , x (!) 

p(r) = Ci n pi mn ji(k n r)Y lm (d,(j>) 



l — m— — l n — 1 



^ Cln pimnji(k n r)Yi m (6, 



(4) 



Here, pi mn is the density coefficient, Yi m (9, (f>) is the spherical 
harmonic corresponding to the spherical angular coordinates, ji (x) 
is the spherical Bessel function and Cj n is the normalisation con- 
stant of the spherical Bessel Function. The resolution of the angu- 
lar modes is determined by the value of Z max and that of the radial 
modes by the value of n max (/). In theory, Z max and n max (Z) tend to 
infinity. In practice, since the distances which galaxy surveys probe 
are finite, Z max and n max (Z) need to be truncated at finite numbers. 
This truncation is an effective smoothing of the density field and 
limits the resolution of the structure in the reconstruction such that 
the structures smaller than the effective resolution of the highest 
harmonic, A8 ~ 7r/Z max (cf. Peebles 1980, Section 46) and the 
highest radial mode, Ar ~ i~/fcn max , will be lost. On the other 
hand, if the expansion is continued to a very high order the resolu- 
tion of expansion will be higher than the size of the real structure, 
introducing artifacts into the reconstruction of density field. The 
choice of Z max and n max (Z) is discussed in Section 6. 

The radial modes, k n , correspond to a given angular mode 
I and the C*; n s are chosen according to the boundary conditions 
and to ensure the orthogonality of the spherical Bessel functions. 
It is assumed that the data are only given within a sphere of ra- 
dius R and that the logarithmic derivative of the gravitational po- 
tential is continuous at r = R. The corresponding value of C; n 
is 1/ (^^j-[ji(k n R)] 2 ^ . The derivation of this result is given in 
FLHLZ, Appendix A. 

Given that the data from a flux-limited redshift catalogue are 
a set of N discrete points with a redshift-space position vector s 
and a weighting function iu(si), the harmonic density coefficients 
in redshift space can be derived using the inversion formula, 
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s 

Plmn 



w(si)ji(k n Si)Y* m (9i, 4>i) 



(5) 



where the superscripts S and R hereafter denote redshift and real- 
space, respectively. We give each galaxy the weight w(si) = 
l/(f) s (si) and 4> s (si)(tt 4> s (ri)) is the selection function of galaxy 
i. 

The density fluctuations, S(s) = (p(s)/p) — 1, can be esti- 
mated in terms of the Fourier-Bessel function 



cS Plmn 

Vlmn — 



M, 



(6) 



where p is the mean density of the survey and the second term 
M; mn represents a monopole correction (cf. FLHLZ, Equation 
D14). 

3.2 Expansion of the velocity field 

The linear theory velocity field in spherical harmonics is given by 
Regos & Szalay (1989). It is derived by expanding the gravita- 
tional potential <f>(r) and the density fluctuations <5(r) in Fourier- 
Bessel series, and substituting these into Poisson's equation. Since 
ji(k„r)Yi m (9, (j>) is an eigenfunction of the Laplacian operator, 



V 2 [ji(k n r)Y lm (e, 



-k n ji{k n r)Y lm {6,4>), 



(7) 



the harmonics of the gravity field can be related to those of the 
galaxy density field as: 



3 (-.0.6 tj2 filmn 



(8) 



The radial component of the velocity field, vy (r) = f • v(r), 
can be decomposed into harmonics as in Equation 4. The harmonic 
coefficients v lmn for the predicted radial velocity field within r < 
R are derived by comparing the harmonic expansions of the veloc- 
ity field and the gravitational potential field: 



Vlr, 



= I3H Y^ (Si)w (*w)> 



(9) 



with H as the velocity matrix (cf. FLHLZ, Equation 22). 

In spherical harmonics, the LG velocity, v(0), is written as 
(see FLHLZ, Appendix CI): 



v(0) 



PH 



3V4tt 



E Ci„ / dr'ji(k n r') 



x (- V2Re[5? ln ] i + >/2Im[6? ln ] j + Re[5? 0n ] k)l0) 

where 5ft [a] and Q[a] refer to the real and imaginary parts of a com- 
plex number, a, and R is the radius out to which the density field is 
measured for the calculation of the velocity. 

The harmonics of the transverse velocity field (v_i_ = — r x 
[r x v(r)]) are also related to those of the density field. They can 
be derived using quantum mechanical methods (see FLHLZ, Ap- 
pendix C2). 



4 DECONVOLVING REDSHIFT-SPACE DISTORTIONS 
IN SPHERICAL HARMONICS 

The galaxy peculiar velocities will only introduce distortions in the 
radial direction, leaving the angular parts intact. Thus, if the density 
field in redshift-space is expressed in terms of the orthogonal radial 
and the angular components, as it is the case in the Fourier-Bessel 



space, the redshift distortions will only couple the radial modes of 
the real space density field. This coupling is described by a cou- 
pling matrix, (Z;) nn , (Fisher et al. 1994 and FLHLZ): 



*(0 



$lmn — E (^)nn' 



(11) 



where (Z;) nn / depends only on Ho, (3, w(r) and <f> B (r) and is in- 
dependent of the power spectrum 5 . The full derivation of the cou- 
pling matrix is given in FLHLZ, Appendix D. We note that Equa- 
tion 1 1 is calculated for linear scales only. On small scales, the pe- 
culiar velocities introduce radially extended distortions called Fin- 
gers of God. These distortions can be corrected for either by col- 
lapsing the fingers to a point at high density regions or by assum- 
ing a Maxwellian distribution for the peculiar velocities (Peacock 
& Dodds 1994) and deconvolving this distribution from the recon- 
structed fields (Heavens & Taylor 1995; Erdogdu et al. 2004). In 
any case, they will not change the constructed fields substantially 
as they are smoothed out due to the resolution of the spherical har- 
monics (see Erdogdu et al. 2004 who use who a similar resolution). 

In the absence of shot-noise, the real-space density harmonics 
can simply be obtained by inverting the coupling matrix in Equa- 
tion 11. However, in the presence of shot-noise, a straight inversion 
becomes unstable leading to inaccurate estimation of the real-space 
harmonics. This problem can be overcome by Wiener filtering. 



5 WIENER RECONSTRUCTION 

In this section, we give the formulae for the Wiener filter method. 
For a full derivation, we refer the reader to Zaroubi et al. (1995) 
and Erdogdu et al (2004). 



5.1 Model Signal and Noise Matrices 

For this analysis, the Wiener-filter reconstructed real-space density 
field is given by: 

(t)wF = E ( S < + N '] _1 )n„' (^L^ *W' -(12) 



For the case w(r) = l/4> s (r), the signal and noise matrices sim- 
plify to 

(S;)„„' ^ P{k n )C^5* n , and 



If 2 1 

(N ( )w = - drr —^-jiiknrjjiik^r) , 
PJ fair) 



(13) 



(14) 



respectively. 

The expected scatter in the reconstructed density, (AS(r)) 2 is 
given by 



(15) 



(A5(r)) 2 = i-^(2/ + l)[(I-F ( )S i ]w 

Inn' 

x Ci n Ci n 'ji(k n r)ji(k n >r) . 

The scatter in the reconstructed radial velocity (Av r (r)) fields 
can also be formulated in Fourier-Bessel space in an analogous way 

5 The coupling matrix has a more complicated dependence on / and m for 
catalogues with incomplete sky coverage. Zaroubi & Hoffman (1996) de- 
rived a similar expression for the coupling matrix in Cartesian coordinates. 
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(cf. FLHLZ, Appendix F): 

(Av r (r)) 2 = ^£(2i + l)[(I-F,)S,W 

Inn' 

r , n j'i{k n r)j'i(k n ir) 

X O in C in / — — . (10) 

fcnfcn' 

5.2 Choice of Reference Frame 

The redshift distortion matrix given in Section 4 takes into account 
the motion of the observer, v(0), and thus the choice of reference 
frame is arbitrary. On the other hand, the formulation of (Z;) nn , 
involves the Taylor expansion of the velocity field out to first order 
in Av = v(r) — v(0) (see Appendix D of FLHLZ). As such, the 
reconstruction will be more accurate in the frame in which Av is 
smaller. Nearby where the galaxies share the LG motion, it is more 
accurate to correct for the redshift distortions in the LG frame. At 
larger distances, galaxy motions are independent of the LG velocity 
and thus Av is smaller in the CMB frame. Since 2MRS probes the 
local universe (cz mc d ~ 6000 kms -1 ), we will work with the LG 
frame redshifts. We will present the dipole velocity results from 
redshifts in both frames. 

The choice of frame for the redshifts will not affect the frame 
of the reconstructed density field. In both cases, the reconstructed 
harmonics will be in real-space. Consequently, the reconstructed 
velocity field (Equation 9) will be in the CMB frame. It is possi- 
ble to convert the velocity field into the LG frame by subtracting 
the Local Group velocity (vlg = 627 ± 22 kms -1 , towards 
Ilg = 273° ± 3°, b LG = 29° ± 3°, Bennett et al. 2003 and 
Courteau & Van Den Bergh 1999). In the following section, the re- 
constructed velocity maps will be presented in either the LG or the 
CMB frames. Up to 8000 kms -1 (Figures 11, 12, 13 and 14) it is 
difficult to see the real nature of the structures in the LG frame as 
the LG velocity dominates. Therefore, the maps will be plotted in 
the CMB frame. On larger distances (Figures 15, 16, 17 and 18), 
galaxies have positive CMB radial velocities so it is easier to com- 
pare the velocities in the LG frame than in the CMB frame. 



6 APPLICATION TO THE TWO MASS REDSHIFT 
SURVEY 

The formalism discussed in the previous sections is limited to 47T 
sky coverage. In principle, the analysis can be extended to account 
explicitly for incomplete sky coverage. In this paper, a simpler 
approach has been adopted to fill in the galaxies masked by the 
Galactic Plane (the Zone of Avoidance). The masked area is divided 
into 36 bins of 10° in longitude. In each angular bin, the distance 
is divided into bins of 10 hT 1 Mpc. The galaxies in each longi- 
tude/distance bin are then sampled from the corresponding longi- 
tude/distance bins in the adjacent strips — b mas kcd — 10° < b < 
l&maskcd! + 10° (where |& mas ked| = 5° or & maskc d = 10°). These 
galaxies are then placed in random latitudes within the mask region. 
The number of galaxies in each masked bin is set to a random Pois- 
son deviate whose mean equals the mean number of galaxies in the 
adjacent unmasked strips. This procedure is carried out to mimic 
the shot noise effects. The success of this interpolation method de- 
pends on the interplay between the width of the mask, the angular 
resolution and whether the structure in the adjacent regions phys- 
ically correlate with the structure in the Zone of Avoidance. The 
method is less robust for masked areas larger than |b| = 15° (La- 
hav et al. 1994). The width of the 2MRS mask is much smaller than 



this value and the survey penetrates deep into the Zone of Avoid- 
ance apart from very obscured regions near the centre of the Milky 
Way. We also test the effects of our method by varying the bin 
sizes and by filling the region uniformly instead of interpolating. 
We calculate the angular power-spectrum for each case and find 
good agreement (Rassat et al. 2006) between them. 

The density field is expanded within a spherical volume of ra- 
dius i? max = 20000 km s _1 . We impose the condition that the 
density field outside i? max is zero and the logarithmic derivative 
of the gravitational potential is continuous at the boundary (which 
holds if ji-i (k n R) — for all I). There are two other possible 
boundary conditions discussed in FLHLZ: setting only the density 
or the radial velocity zero at the boundary. We have tested our re- 
constructions using these boundary conditions. Setting the density 
to zero at the boundary results in a discontinuous potential and the 
reconstructed radial velocities become unrealistically high. Up to a 
distance of 10000 kms -1 , the fields reconstructed using the zero 
velocity boundary condition agree reasonably well with those re- 
constructed using the continuous potential condition, however the 
radial velocities of the zero velocity reconstructions are slightly 
higher. Beyond 10000 kms -1 , the radial velocities of the zero 
velocity reconstructions tend to zero, however the backside infall 
into the Shapley Supercluster still exists. Thus we conclude that 
setting the logarithmic derivative of the gravitational potential to 
be continuous at the boundary allows for a smooth transition of the 
density and velocity fields from r < i? max to r > i? max and as 
such gives the most robust results. We also note that FLHLZ, WLF 
and Schmoldt et al. (1999) use the same boundary condition for 
their reconstructions. 

The resolution of the reconstruction depends on the maximum 
values of angular and radial modes (see Figure 6 of WLF). We want 
to have as high a resolution as possible without introducing artifi- 
cial structures to the field. Therefore, we compute the harmonics 
up to £ max = 15 angularly and k„R ^ 100 radially. Z max = 15 is 
chosen so that the angular resolution (59 ~ n/l) is comparable to 
the width of the interpolated region and the radial modes are cho- 
sen to match the angular resolution (see FLHLZ, Appendix B for a 
detailed discussion). 

The Wiener filter method requires a model for the linear 
galaxy power spectrum in redshift-space, which depends on the 
real-space power spectrum and on the redshift distortion param- 
eter, (3 = n m °' 6 /6. The real-space power spectrum is well de- 
scribed by a scale-invariant Cold Dark Matter power spectrum with 
shape parameter Y w Q. m h, for the scales concerned in this anal- 
ysis. The combined analysis of third year data from the Wilkinson 
Microwave Anisotropy Probe (WMAP, Spergel et al. 2006), the 
Two Degree Field Galaxy Redshift Survey (2dFGRS, Cole et al. 
2005) and the Sloan Digital Sky Survey (Tegmark et al. 2004) im- 
ply T « 0.17. On the other hand Frith, Outram & Shanks (2005) 
and Mailer et al. (2005) analyse the angular power spectrum of the 
2MASS galaxies and derive lower F values of 0.14 ± 0.02 and 
0.116 ± 0.009, respectively. The normalisation of the power spec- 
trum is conventionally expressed in terms of the variance of the 
mass density field in spheres of 8/iT 1 Mpc, as. Recent WMAP 
results combined with the 2dFGRS gives o$ w 0.75, whereas anal- 
ysis of the 2MASS galaxies suggest ag « 0.9 (Frith, Outram & 
Shanks 2005, Mailer et al. 2005, Pike & Hudson 2005). In a pre- 
vious paper, we analysed the 2MRS LG dipole and derived a value 
for /3 for 2MRS: f3 = 0.4 ± 0.09 (Erdogdu et al. 2006). 

In order to check the goodness of fit of these different param- 
eters to the data, we calculate the \ 2 statistic defined by 
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X 2 =d t (S + N)- 1 d, (17) 

for T = (0.1,0.2), a 8 = (0.7, 0.9) and /3 = (0.3,0.4,0.5). 
Changing V does not affect the \ 2 values greatly but V = 0.2 
is a slightly better fit. Varying as is a bit more influential on the 
X 2 values, however, the reconstructed maps look very similar for 
these models. Increasing F decreases the power on large scales and, 
therefore, the peculiar velocities tend to be smaller but the differ- 
ence is negligible. Decreasing as has the same effect on the power 
spectrum as increasing T. The choice of /3 plays a more important 
role in the goodness-of-fit results but has little affect on the recon- 
structed maps other than a linear scaling of peculiar velocities. In 
the analysis that follows, we adopt T — 0.2, as = 0.7 and f3 = 0.5, 
which is the best fit model that was considered, with a reduced \ 2 
value of 1.3. We note that the best fit f3 is higher than derived in 
Erdogdu et al. (2006), this discrepancy will be discussed in Section 
6.3. In a forthcoming paper (Rassat et al. 2006), we will present a 
detailed analysis of the cosmological parameters of 2MRS. 



6.1 Density Maps 

Figures 3, 4, 5, 6, 7, 8, 9 and 10 show the Aitoff projections of the 
reconstructed density field in real-space, plotted in Galactic coor- 
dinates, evaluated at different distance slices. 

Figure 3 shows a thin density shell at r — 2000 km s _1 . The 
expected scatter in the reconstructed density field is (A<5(r)} = 
0.31 and this value stays approximately the same out to r = 
10000 km s -1 and increases by 0. 1 thereafter. All the known struc- 
tures are resolved and the field at this distance is dominated by 
inter-connecting voids. The Puppis cluster (Lahav et al. 1993, SI 
in Saunders et al. 1991) is very prominent and extends down to 
b ~ -28°. The overdensity labelled as CI (I ~ 190° ,b ~ 5°) 
does not correspond to any previously observed structure. 

Figure 4 is the 2MRS density field at r = 4000 km s~ 1 . Hy- 
dra (I ~ 270°, b ~ 50°) and Centaurus (/ ~ 300°, b ~ 20°) 
continue in this shell but now are separated by the Hydra void. 
Centaurus is connected to the new overdensities C2 (I ~ 345°, 
b ~ 10°) and C8 (/ ~ 295°, b ~ 7°) and Hydra connects 
to C4 (I ~ 240°, b ~ 25°). Pavo-Indus-Telescopium continues 
deep into the Zone of Avoidance and connects to the Centaurus 
Wall as suggested by Fairall (1988). The well-known cluster Can- 
cer (I ~ 210°, b ~ 30°), as well as new clusters C3 (I ~ 275°, 
b ~ -30°), C5 (I ~ 195°, b ~ 0°), C6 (/ ~ 170°, b ~ -15°) 
and C7 (I ~ 55°, b ~ 15°), all peak at r = 5000 km s" 1 with C5, 
C6 and C7 connecting to Pegasus and the Perseus-Pisces Wall. The 
overdensity Ce (WLF) is more prominent in the 2MRS field than in 
IRAS 1.2 Jy field. The Perseus and Pisces clusters are picked up 
as two distinct overdensities and both connect to Camelopardalis 
(/ ~ 145° ,6 ~ 30°) and Cetus (I ~ 190°, b ~ -55°) at 
r = 5000 km s -1 . 

Figure 5 is the density field at r = 6000 km s -1 . Hydra 
(I ~ 285°, b ~ 5°) is less prominent in this shell whereas Pavo- 
Indus-Telescopium (I ~ 0°, b ~ —45°) is still strong. The cluster 
Abell 3627 (or Norma, I ~ 330°, b ~ -10°) which was recog- 
nised as the centre of the Great Attractor 6 by Kraan-Korteweg et 
al. (1996), peaks around r = 5000 km s -1 and is separated from 
the Centaurus Wall by the Sculptor Void. The overdensity labelled 



By 'Great Attractor', it is meant the entire steradian on the sky centred at 
(I ~ 310° ,b ~ 20°) covering a distance of 20 h' 1 Mpc to 60 h' 1 Mpc. 



as C8 continues from the previous shell and contains the lumi- 
nous CIZA J1324.7-5736, cluster which was discovered by Ebel- 
ing, Mullis & Tully (2002) and the Centaurus-Crux cluster (Woudt, 
1998). The Abell clusters 569 (Z ~ 155°, b ~ -15°) and 168 
(I ~ 120°, b ~ 70°) are as prominent as the Perseus-Pisces ridge 
and both are connected to it. The clusters labelled as WLB (e.g. 
WLB 248 at I ~ 240°, b ~ 40°) were identified by White et al. 
(1999). There are two previously unidentified voids at this distance 
labelled as VI and V2 and a cluster at (Z ~ 240°, b ~ 8°) labelled 
C9. 

Figure 6 is the 2MRS density field at r = 8000 km s -1 . 
The Great Wall (Geller & Huchra 1989) connects to the Ophiuchus 
supercluster discovered by Wakamatsu et al. (1994). The massive 
Shapley supercluster is already visible (I ~ 300°, b ~ 30°) at this 
distance. The void VI continues to this shell. The new clusters C10 
(I ~ 205°, b ~ 10°), Cll (Z ~ 250°, b ~ -10°) , C12 (Z ~ 325°, 
b ~ -10°) and C13 (/ ~ 85°, b ~ 13°) lie on galactic latitudes. 
The cluster ZwCl 0943.7+5454 was first identified by Zwicky & 
Herzog (1966) and K42 by Klemola (1969). All the superclusters 
labelled as SC (e.g. SC 129 in this shell) are given in Einasto et al. 
(1997, Table 2) and structures labelled CAN and CID (e.g. CAN 
136 in this shell and CID 15 in the next shell) are clusters given in 
WegnereraZ. (1999). 

Figure 7 is the 2MRS density field at r = 10000 km s _1 . 
The expected scatter now rises to (AS(r)) = 0.36. The mighty 
Hercules supercluster dominates this shell and Shapley begins to 
appear. Other dominant structures are CID (CAN) 15 & 16. 

Figure 8 is the 2MRS density field at r = 12000 kms -1 and 
has a scatter of (A<5(r)} = 0.41. The Perseus-Pegasus and Pisces 
(Einasto etal. 1997) superclusters, C20 and Abell 576 continue the 
Perseus-Pisces wall and dominate the shell. 

Figure 9 is the 2MRS density field at r = 14000 km s _1 
with scatter (A<5(r)} = 0.46. The density field of the Shapley su- 
percluster peaks at this distance and connects to C28 and C29 at 
r = 15000 km s~ 1 to form a ridge. 

Figure 10 is the 2MRS density field at r = 16000 kms -1 . 
The expected scatter is (A<5(r)) = 0.49. The most conspicuous 
structure is the Pisces-Cetus supercluster, which connects to the 
Southern Great Wall. The Perseus-Cetus also dominates the south- 
ern strips of the 2dF (see e.g. Erdogdu et al. 2004, Porter & Ray- 
chaudhury 2005) and Sloan surveys (e.g. Porter & Raychaudhury 
2005). We note that at this distance, we do not pick up any voids. 
This is due to the nature of our technique. The Wiener filter method 
will predict the mean field in the absence of data. Since the den- 
sity field is constructed to have zero mean, the Wiener filter signal 
will approach to zero at the edges of the survey where the shot 
noise dominates. This leads to the signal being reconstructed in a 
non-uniform manner. As such, we do not present maps of the den- 
sity and velocity fields beyond 16000 kms -1 after which the shot 
noise is too high and the Wiener Filtered fields tend to zero. 

6.2 Velocity Maps 

Figures 11, 12, 13, 14, 15, 16, 17 and 18 show the Galactic 
Aitoff projections of the radial velocity fields for /3 = 0.5 eval- 
uated across the same shells as the density fields shown in Fig- 
ures 3, 4, 5, 6, 7, 8, 9 and 10. Positive (outflowing) radial ve- 
locities are shown as solid lines and the negative (Mailing) ra- 
dial velocities are shown as dashed lines. The contour spacing is 

\A-Vradial\ = 50kmS _1 . 

Figures 11 & 12 show the CMB frame radial velocity fields at 
r = 2000 kms -1 and at r = 4000 km s~ 1 , respectively. The ex- 
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Figure 3. The reconstructed density field, evaluated on a thin shell at 2000 km s , shown in Galactic Aitoff projection. Dashed lines show S < 0, and 
solid lines show 8^0, with contour spacing |A<5| = 0.1. The overdense regions are (from left to right) Ursa Major (Urs-Maj),WLB 550, Co, Virgo (Vir), 
Centaurus (Cen), Hydra (Hyd), Puppis (Pup), Fornax-Doradus-Eridanus (F-D-E) and CI . The voids are Gemini (Gem), Taurus (Tau), Andromeda (And), 
Delphinus (Del), Virgo (Vir), Local Void (LV), Eridanus (Erid) and Orion (Or). 



Bfjf 
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Figure 4. Same as Figure 3 but evaluated at 4000 km s 1 . The overdense regions are C6, Perseus (Per), Pisces (Pis), Camelopardalis (Cam), C/3, Pegasus 
(Peg), CS, Cj, Ce, C7, Cf, C2, Coma (Com), Virgo (Vir), Centaurus (Cen), Pavo-Indus-Telescopium (P-I-T), C8, Hydra (Hyd), C3, C4, Cancer (Can), NGC 
1600 (N1600), Cetus (Get) and C5. The voids are Aunga (Aun), Taurus (Tau), Perseus-Pisces (Per-Pis), Cygnus (Cyg), Aquarius (Aqu), Corona Borealis 
(Cor-Bor), Microscopium (Mic), Virgo (Vir), Sagittarius (Sag), Hydra (Hyd), Octans (Oct), Leo, Eridanus (Erid) and Canis Major (Can-Maj). 
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Figure 5. Same as in Figure 3 but evaluated at 6000 km s -1 . The overdensities are Abell 569 (A569), Perseus (Per), Pisces (Pis), NGC 1600 (N1600), Abell 
168, Camelopardalis (Cam), Pegasus (Peg), C7, C£, Pavo-Indus-Telescopium (P-I-T), Abell Clusters 1605, 300 & 42, Abell 3627, C8, Phoenix (Pho), Abell 
2923, Coma (Com), Hydra (Hyd), WLB 248, C9, Cancer (Can), Orion (Or) and Cetus (Cet). The voids are Taurus (Tau), Ursa Major (Urs-Maj), Perseus-Pisces 
(Per-Pis), Aquarius (Aqu), Corona-Borealis (Cor-Bor), Pegasus (Peg), V2, Microscopium (Mic), Coma (Com), Sculptor (Scu), Columba (Col), VI and Lepus 
(Lep). 




Figure 6. Same as in Figure 3 but evaluated at 8000 km s -1 . The overdensities are ZwCl 0943.7+5454 (ZWI), Abell 634, Perseus-Pisces (Per-Pis), Abell 
400, WLB 664, C13, Pegasus (Peg), Abell Clusters 2634 & 2666, Abell S0929, Klemola 42 (K42), Great Wall (GW), Ophiuchus (Oph), SC 129, C12, Phoenix 
(Pho), Coma (Com), SDSS CE J159.778641-00.784376 (SDSS), Abell 3389, WLB 213, Cll, CAN 136, C10 and Abell 539. The voids are Perseus-Pisces 
(Per-Pis), Fornax (For), Serpens (Ser), Capricomus (Cap), Apus (Ap), Columba (Col) and VI. 
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Figure 7. Same as in Figure 3 but evaluated at 10000 kms -1 . The overdensities are CID 15, CID16, C13, Abell 2634, C14, CAN 75, CAN 338, Hercules 
(Her), SC180, Shapley (Shap), Phoenix (Pho), Leo, RXC J07 12.0-6030, Lepus (Lep), ZwCl 0820.6+0436, Abell 496 and CIZA J0603.8+2939. The voids are 
V9, V6, V10, V12, VI 1, V13, V5, V4, V3, Apus (Ap), V7, Columba(Col) and V8. 
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Figure 8. Same as in Figure 3 but evaluated at 12000 kms -1 . The overdensities are SDSS clusters SDSS CE J048.551922-00.613762 and SDSS CE 
J050.935932-00. 117944 (SDSS), Abell 576, C20, Pisces (Pis), CAN 137, C19, C18, Perseus-Pegasus (Per-Peg), C16, C15, Hercules (Her), C17, SC180, 
SC194, Shapley (Shap), C24, C23, Leo, SC73, Abell S0463, C22, Lepus (Lep) and C21. The voids are V6, V14, V5, V4, Further-Fornax (Fur-For), V3 and 
Columba (Col). 
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Figure 9. Same as in Figure 3 but evaluated at 14000 kms . The overdensities are Rixos F231.526 (RIX), SC44, C19, Pisces (Pis), Perseus-Pegasus 
(Per-Peg), C25, C26, Hercules (Her), Abell S0757, C28, Shapley (Shap), C29, C30, SC 43, Leo, Abell 3376, C27 and C21. The voids are V4, Further-Fornax 
(Fur-For)andV15. 




Figure 10. Same as in Figure 3 but evaluated at 16000 kms -1 . The overdensities are Abell 582, C31, Ursa-Major (Urs-Maj), SC168, West 9 (W9), C33, 
Abell 2319 & Cygnus A (CYGA), Pisces-Cetus (Pis-Cet), C32, Aquarius-Cetus (Aq-Cet), Hercules (Her), SC143, Southern Great Wall (SGW), C29, Shapley 
(Shap), C30, SC43, Sextans (Sex) and SC68. 
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pected scatter in the velocity field is ( Av(r)} = 53 km s 1 at r — 
2000 kms" 1 and (Av(r)> = 62kms" 1 at r = 4000 kins" 1 . 
At r = 2000 kms" 1 , there is a strong outflow towards the Virgo- 
Great Attractor (I ~ 290°, b ~ 15°) and inflow out of the 
underdense regions. Also seen at this distances are the effects 
of Perseus-Pisces and N1600, causing a strong outflow towards 
(I ~ 165°, b ~ -30°). At r = 4000 kms" 1 , Perseus-Pisces 
and N1600 are the dominant structures causing an outflow towards 
(I ~ 150°, b ~ —20°). There is also a flow towards the centre 
of Pavo-Indus-Telescopium (Z ~ 335°, b ~ —25°) and a back- 
side infall to Hydra. The back-side infall to the Great Attractor 
region becomes strongest at r = 5000 kms" 1 and reaches to 
finfaii = (982 ± 400)/3kms" 1 in the Local Group frame and 
Winfaii = (127 ± 409)/° kms" 1 in the CMB frame, where the 
quoted values are the mean of the velocity field around the region 
and the error bars are the standard deviation from that mean and the 
expected scatter in the velocity field added in quadrature. 

Although still visible at r = 4000 kms" 1 , the Great 
Wall really starts to dominate in Figure 13 at a distance of 
r = 6000 kms" 1 . At this distance the velocity scatter rises to 
(Av(r)} = 72 kms -1 . Here the fall towards the Coma Clus- 
ter is quite apparent at I ~ 240° and b ~ 80°. The velocity 
flow towards Perseus-Pisces is still strong and there is some out- 
flow towards Pegasus (7 ~ 95°, b ~ -60°), Phoenix (I ~ 330°, 
b ~ -75°), the Abell clusters A548, A539 (I ~ 230°, b ~ -20°) 
and A400 (I ~ 175°, b ~ -45°). The CMB radial velocities at 
r — 8000 km s" 1 , shown in Figure 14, have scatter of (Av(r)} = 
87 kms" 1 and are mostly outflowing apart from the remnants of 
back-side infall towards the Great Attractor and the Coma clus- 
ter. By r = 10000 kms -1 , it becomes difficult to compare pecu- 
liar velocities in the CMB frame so we switch to the LG frame. 
Atr = 10000 kms -1 (Figure 15, (Av(r)) = 108 kms -1 ), the 
dipole velocity of the LG is still visible, however there is an outflow 
towards Shapley (I ~ 310°, b ~ 25°) and the Hercules superclus- 
ter (I ~ 20°, b ~ 45°). Shapley outflow becomes stronger at r = 
12000 kms -1 (Figure 16, (Av(r)) = 130 km s _1 ), begins to de- 
crease atr = 14000 kms -1 (Figure 17, (Av(r)) = 146 kms -1 ) 
and atr = 16000 kms" 1 (Figure 18, (Av(r)) = 150 kms" 1 ) 
there is a strong back-side infall towards it. 

Figure 19 shows the 2MRS density and velocity fields in 
the Supergalactic plane (SGZ — Okms" 1 ) in different reference 
frames. The top plot is the reconstructed three-dimensional pecu- 
liar velocity field in the CMB frame, superimposed on the recon- 
structed density field in the Supergalactic plane. The bottom plot 
is the same as the one on the left but the superimposed peculiar 
velocity field is reconstructed in the LG frame (obtained by sub- 
tracting the Local Group velocity (vlg = 627 ± 22 kms" 1 , to- 
wards I lg = 273° ± 3°, b LG = 29° ± 3°). The extent of the 
planar density distribution in the Supergalactic plane can clearly 
be seen in these plots. Since the Local Group is moving towards 
the Great Attractor/Shapley region, the strong outflow towards the 
positive y-plane seen in the top plot decreases in the bottom plot. 
For the same reason, the outflow towards the negative y-plane is 
much stronger in the LG frame (bottom plot) than in the CMB 
frame. The most intriguing result is the connectivity of the Shap- 
ley supercluster (SGX,SGY) = (-12000, 7000) to the Great At- 
tractor region. The Virgo cluster at (SGX, SGY) = (0, 1000) 
also appears connected to the Great Attractor. The Perseus-Pisces 
supercluster also is visible at (SGX,SGY) = (2500,-4500) and 
Camelopardalis is at (SGX,SGY) = (45, -30). The other visible su- 
perclusters are Columba at (SGX,SGY) = (7500, -11000) and Her- 
cules (SGX,SGY)=(- 13000, -4000). In both plots, the outflow is 



dominated by Shapley, particularly in the CMB frame. The flow to- 
wards Virgo and the Perseus-Pisces region is also apparent and the 
flow towards the Great Attractor is clearly visible. Also observed 
is a significant flow towards Columba and some towards Hercules. 
After 16000 kms -1 , the shot noise is too high and the density 
field becomes non-uniform. 

Over the past years, there has been some disagreement over 
the existence of back-side infall towards the GA (e.g. Dressier & 
Faber 1990 and Mathewson et al. 1992). We find a clear back-side 
infall, visible in both plots and in particular the one in the LG frame. 
The highest velocity back-falling velocity around the GA region is 
800 ± 600/T 1 kms" 1 (at 55ft" 1 Mpc) in the CMB frame. In 
Figure 20, we plot the average radial velocity field in towards the 
centres of the GA (290° s? / ^ 320°, -25° < b < 60°-25° < 
b sC 70°, r < 80ft" 1 Mpc ) and the Shapley (290° < I < 320°, 
40° < b < 70°, r > 80ft" 1 Mpc) Superclusters as a function 
of distance. The red lines are for the LG frame and the black lines 
are for the CMB frame. The solid lines denote the mean of the 
reconstructed velocities and the dashed lines are the errors calcu- 
lated from the standard deviation of velocities in that region and 
the expected velocity scatter added in quadrature. In this plot, the 
back infall onto the GA and the subsequent outflow towards Shap- 
ley are clearly visible. An observer in the GA would see the sur- 
rounding nearby galaxies moving towards her. On larger scales, she 
would observe a dipole velocity towards Shapley. If she was in the 
Shapley Supercluster (at r « 130ft" 1 Mpc), she would observe 
that the surrounding galaxies are moving towards her. Whilst, this 
back-side infall onto the GA was also observed by FLHLZ, WLF 
and Schmoldt et al. (1999), other authors who applied different re- 
construction techniques to IRAS galaxy catalogues (e.g. Davis, 
Nusser & Willick 1996; Branchini et al. 1999; Valentine, Saun- 
ders and Taylor 2000) find little evidence of it. Although the 2MRS 
sample is the densest all-sky sample to date and as such allows the 
most precise reconstruction of the local velocity field, our results 
are still not conclusive. In a forthcoming paper we will compare 
the reconstructed fields with the observed velocity fields based on 
the distance measurements to elliptical galaxies in the 6dFGS and 
a new survey of spirals in the I Band (SFI++) to assess the signifi- 
cance of this back-side infall. 

6.3 Alternative Reconstruction Methods 

In order to assess the effects of the systematic errors of our method, 
we also carried out the reconstructions using alternative methods. 
Figure 21 shows several of these reconstructions on the Super- 
galactic Plane (SGP) out to 6000 kms" 1 . The top left plot is the 
unsmoothened harmonic reconstruction of the 2MRS density field 
in redshift space. It can clearly be seen that the harmonics are 
strongly affected by shot noise. Overdensities follow circular paths 
about the origin as a result of correlations introduced by noise. 
These artifacts are still clearly evident in the top right plot where 
the density field is corrected for redshift distortions (using Equa- 
tion 11) but without noise supression. The correction of the red- 
shift distortions reduces the overall amplitude of the density con- 
trasts and makes the structures more spherical. This reconstruction 
corresponds to the unbiased minimal variance (umv) estimation of 
Zaroubi (2002) since in our case the redshift distortion matrix is 
invertible. The bottom left panel shows the real space density field 
after the Wiener Filter has been applied. As expected the Wiener 
Filter mitigates the effects of shot noise, removing the artifacts and 
high frequency harmonics from the map, however, the Wiener re- 
constructed field is biased towards the mean field which is zero by 
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Figure 11. The reconstructed radial velocity field, evaluated on a thin shell at 2000 km s~ 1 , shown in Galactic Aitoff projection in the CMB frame. This 
velocity field corresponds with the density field shown in Figure 3. Dashed lines show infall, and solid lines show outflow. The first solid line is for ^ ra dial = 
Okms -1 , and contour spacing is |Av radial | = 50kms _1 . Note the strong outflow towards the Great Attractor and flow out of the Local Void. Also 
Perseus-Pisces and N1600, cause an outflow towards (I ~ f65°, b ~ —30°) 




Figure 12. Same as Figure 1 1 but evaluated at 4000 km s~ 1 . This velocity field corresponds with the density field shown in Figure 4. Note Perseus-Pisces and 
N1600 cause the strongest outflow towards (I ~ f 50°, 6 ~ —20°). There is also a flow towards the centre of Pavo-Indus-Telescopium (I ~ 335°, b ~ —25°) 
and a weaker one towards Hydra (I ~ 280°, b ~ 10°). 
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Figure 13. Same as Figure 1 1 but evaluated at 6000 kms -1 . This velocity field corresponds with the density field shown in Figure 5. Note the strong outflow 
towards the Great Wall, especially towards Coma and the general velocity flow towards the North Galactic Pole from the South Galactic Pole. There is still 
some outflow towards Pegasus. Also visible is a strong back-side infall towards the Great Attractor region and Perseus-Pisces. 




Figure 14. Same as Figure 11 but evaluated at 8000 km s -1 . This velocity field corresponds with the density field shown in Figure 5. There is still a strong 
outflow towards the clusters of the Great Wall. The outflow towards Shapley (I ~ 310°, 6 ~ 25°) is already visible. The Sculptor Wall is implied in the 
outflow towards the Southern Galactic Pole. 
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Figure 15. Same as Figure 1 1 but evaluated at 10000 km s 1 and in the LG Frame instead of the CMB Frame. There are three superclusters that dominate 
the outflow, Hercules towards (I ~ 20°, b ~ 45°), Shapley towards (/ ~ 320°, b ~ 25°) and Columba towards (/ ~ 225°, b ~ -30°). 
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Figure 16. Same as Figure 1 1 but evaluated at 12000 km s 1 and in the LG Frame instead of the CMB Frame. 
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Figure 17. Same as Figure 1 1 but evaluated at 14000 km s 1 and in the LG Frame instead of the CMB Frame. 




Figure 18. Same as Figure 1 1 but evaluated at 16000 km s 1 and in the LG Frame instead of the CMB Frame. 



construction. In fact, the variance of the Wiener Filtered field is al- 
ways smaller than the variance of the true underlying field. In order 
to overcome this drawback, Yahil (1994) introduced an alternative 
filter which preserves the power in the underlying field. This filter 
is defined as the square root of the Wiener Filter. With this defini- 



tion, the filtered field equals to that of the true field. The bottom 
right panel shows the real-space density field reconstructed using 
Yahil's filter. Whilst the estimated field is unbiased, the reconstruc- 
tion is not optimal in the sense of minimum variance and there are 
artifacts due to insufficient reduction of noise. 
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§^^W6()^ftA§??vfSfl ! ?A ! S r B6lf c l)§8BSte8 ^ e '^ s m tne su P er galactic plane. The spacing of the density contours is |A<5 = 0.1| with dashed lines denoting the 
negative and the solid line denoting the positive density contrast. Superimposed in black are the reconstructed three-dimensional peculiar velocity field in 
the CMB (top plot) and the LG (bottom plot) frames. The length of the arrows are about 300 kms" 1 per cell. The main overdensities are Hydra-Centaurus 
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Figure 21. Recontructions of the 2MRS density field in the Supergalactic Plane. All contours are spaced at | A<5| = 0.5, with solid (dashed) contours denoting 
positive (negative) contours. Top left: The density field expanded in spherical harmonics with no smoothing. Top right: Same as in top left but with the 
harmonics corrected for redshift distortion. Bottom left: Same as in top right but smoothed by the Wiener Filter. Bottom Right: The real space density field 
smoothed by the power preseving filter. 



6.4 The acceleration of the Local Group 

A detailed calculation of the acceleration of LG using the 2MRS 
data was given by Erdogdu et al. (2006). In this section, the LG 
velocity is reconstructed from the density field. There are two main 
advantages of the analysis performed in this paper. Firstly, working 
with spherical harmonic allows for reliable correction of the ve- 
locity flows and secondly, the Wiener filtering suppresses the noise 
self-consistently. The LG dipole is linear in the galaxy density field 
(assuming linear biasing). Consequently, the dipole computed from 



the Wiener filter algorithm will be an optimal (in the sense of min- 
imum variance) estimate of the dipole due to matter within the re- 
construction volume as long as non-linear effects can be neglected. 

The reconstructed 2MRS dipole (Equation 10) and the ex- 
pected scatter (Equation 16) in the reconstruction is plotted as a 
function of distance in the top plot of Figure 22 . The bottom plot of 
Figure 22 shows the angle between the direction of the LG dipole 
and the CMB dipole (/ = 273°, b = 29°). The black line is the 
dipole reconstructed with f3 = 0.5 (rescaled to /3 = 1.0) and the 
dashed line denotes the scatter in the velocity reconstruction. The 
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Figure 20. The radial velocity field in the region of the GA and the Shapley 
Superclusters as a function of distance. The red lines are for the LG and 
the black lines are for the CMB frame. The solid lines denote the mean of 
the reconstructed velocities and the dashed lines are the errors calculated 
from the standard deviation of velocities and the expected velocity scatter 
added in quadrature. The regions with negative velocities move towards us 
and positive velocities move away from us. 



inversion of the radial coupling matrix introduces a non-linear re- 
lationship between the value of f3 and the amplitude of the recon- 
structed acceleration. To illustrate this, we also plot the dipole re- 
constructed with ft = 1.0 (the green line). The amplitudes of the 
dipole reconstructions using j3 — 0.5 and j3 = 1.0 are in good 
agreement (within the expected scatter) and the misalignment an- 
gles are almost the same up to 13000 kms -1 . The blue line is 
the dipole reconstructed with f3 = 0.5 from the CMB redshifts. As 
expected, the amplitude of the dipole is underestimated at nearby 
distances and the misalignment angle is bigger than the dipole re- 
constructed from the LG redshifts. We conclude that, since the LG 
dipole is dominated by nearby galaxies, it is more accurate to use 
the LG redshifts in the reconstruction. 

We also plot the number-weighted dipole obtained using raw 
redshifts (solid red line Erdogdu et al. 2006) and the associated shot 
noise (dashed red line). The bottom plot indicates that the direction 
of the Wiener-reconstructed dipole is consistent with the dipole ob- 
tained in Erdogdu et al. 2006). Both dipoles show the 'tug of war' 
between the Great Attractor and Perseus-Pisce and the dipole ve- 
locities are dominated by structure within the distance of 50 h~ 
Mpc. However, the amplitude of the Wiener dipole is less than the 
number-weighted dipole. This is due to the fact that in the previ- 
ous paper, we did not filter the shot noise errors. Furthermore, in 
the previous paper, we did not take peculiar velocities into account 
when calculating the number-weighted dipole and these will lead to 
the over-estimation of the LG-frame dipole (the Rocket Effect). The 
value for /3 we obtain by equating the CMB dipole measurements 
with our reconstructed dipole reflects the decrease in the ampli- 
tude. We get f3 = 0.54 ± 0.12 at 13000 kms -1 whereas we find 
(3 = 0.40 ± 0.09 at 13000 kms" 1 in Erdogdu et al. (2006). The 
convergence of the direction of the LG dipole does not seem af- 
fected by the velocity corrections but there is in general better align- 
ment between the Wiener-reconstructed LG dipole and the CMB 
dipole than the LG dipole, calculated directly from the raw data 
and the CMB dipole. The misalignment angle of the Wiener dipole 



from the LG redshifts is almost as that of the raw-data dipole up to 
13000 km s _1 and the angle is 39° at this distance. At least half of 
this misalignment is probably due to the fact that when we modeled 
the LG dipole, we assumed all the galaxies have the same weight 
(see Erdogdu et al. 2006). Also our analysis uses linear perturbation 
theory which is correct only to first order 8. There may be contribu- 
tions to the LG dipole from small scales which would cause gravity 
and the velocity vectors to misalign. 



7 SUMMARY AND CONCLUSIONS 

In this paper, the Wiener reconstruction technique was applied to 
recover real-space density and velocity fields from 2MRS. The ob- 
served density field was first expanded in terms of spherical har- 
monics and spherical Bessel functions (Equations 4 and 6) which 
are both orthogonal and which together satisfy Poisson's equation. 
Then, the redshift distortions, described by the coupling matrix Z, 
were corrected for by the matrix inversion of Z in Equation 11, 
assuming linear theory and a value of (3 = Q, ' 6 /b = 0.5. Af- 
ter the inversion, the real-space density field was estimated using 
the Wiener filter for a given CDM power spectrum with erg = 0.7 
and r = 0.2 (Equations 12, 13 and 14). Finally, the velocity field 
is reconstructed from the Wiener-filtered density field using Equa- 
tions 9. 

The reconstructed density field resolves most of the known 
structures as well as new clusters and voids. The density and veloc- 
ity fields agree well with the IRAS PSCz and 1.2-Jy fields derived 
using the Wiener filter technique. There is a back-side infall to- 
wards the Great Attractor region, in agreement with Schmoldt et 
al. (1999) and FLHLZ. There is a strong outflow towards the Shap- 
ley supercluster which suggests that Shapley plays an important 
dynamical role in the local velocity field. 

The LG dipole derived from the reconstructed velocity field 
is lower in amplitude at all distances than the LG dipole derived 
from the raw data by Erdogdu et al. (2006). The misalignment an- 
gle between the CMB dipole and the 2MRS dipole decreases down 
to 13° at 5000 kms" 1 and its amplitude is about 800 /3kms _1 . 
If our canonical value for /3 is correct then this suggests most of 
the acceleration is generated within 5000 kms -1 and the structure 
at higher distances contributes less than 30% to the LG dipole. The 
'tug-of-war' between the Great Attractor and the Perseus-Pisces su- 
percluster to which the dipole is commonly attributed can clearly 
be seen in Figure 1 1 . The 2MRS dipole agrees well with the best-fit 
value for the canonical model to the data with j3 — 0.5, for fixed 
T = 0.2 and cr 8 = 0.7. 

The Wiener filter technique performs very well and provides 
a rigourous reconstruction method determined by the data and the 
expected signal. In a forthcoming paper, we will compare the re- 
constructed fields with the peculiar velocities from the 6dF survey. 
This will allow us to test the validity of our reconstructions and 
enable the refinement of our method. 
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Figure 22. Top Plot: The growth of acceleration of the Local Group due to galaxies within a series of successively larger concentric spheres centred on the 
local group. The growth of estimated shot noise is also shown as dashed lines. The black line is the dipole reconstructed with (3 = 0.5. The blue line is the 
dipole reconstructed with /3 = 0.5 from the CMB redshifts. The green line is the dipole reconstructed with /3 = 1.0. The red line is the dipole from Erdogdu 
et al. (2006) and the dashed line is the associated shot noise. Bottom: The angle between the CMB dipole (I = 273°, b = 29°) and the reconstructed Local 
Group velocity. The lines are colour coded to coincide with the top plot. The dashed lines denote the shot noise errors. 
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